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EFFECT OF SUBLIMATION ON THE BEHAVIOR OF CRYODEPOSITS 

by J. Perry Campbell 
Lewis Research Center 

SUMMARY 

The rates of energy loss due to sublimation have been calculated for water, carbon 
dioxide, argon, krypton, xenon, and oxygen over the temperature range in which these 
materials are solids and have vapor pressures of less than 1 torr. These energy losses 
are compared with the possible thermal radiation losses of such deposits in the same 
temperature range. 

The energy loss rate from these cryodeposits due to sublimation is quite sensitive 
to the surface temperature of the deposit. There is only a narrow temperature range 
over which the sublimation and radiation energy losses are of the same magnitude and 
in which both phenomena must be considered in a heat balance. At temperatures above 
this range sublimation is the dominant energy-loss mechanism; at lower temperatures 
the radiation energy loss is the more important. 

The surface temperature, mass loss, and deposit lifetime of a subliming cryo- 
deposit, both with and without an incident energy heat load, are evaluated. Numerical 
results are given for a typical ice cryodeposit, both with and without an incident energy 
heat load from a 300 K blackbody, for substrate temperatures ranging from 190 to 
270 K, and for deposit thicknesses ranging from 10 micrometers to 1 centimeter. 


INTRODUCTION 

Early investigations of cryodeposit properties (refs . 1 to 5) were concerned mainly 
with the basic study of cryopumping and the storage of cryogenic fluids. Later, it was 
recognized that cryodeposits might present a problem to designers of space simulators . 
The cryodeposits formed on the walls of a space simulator, as a result of cryopumping, 
could change the thermal radiative properties of the simulator thereby making it pos- 
sible for radiation to be reflected back to the model. Thus, the totally absorbing heat 
sink of space that must be simulated could be degraded. 

A similar type of argument applies to the storage of cryogenic fluids. A cryodeposit 
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that forms on a cryogenic storage tank, on a shadow shield protecting such a tank, or 
within the multilayer insulation around the tank could alter the thermal radiative proper- 
ties of the system and upset the designed heat balance. 

Actual space flights have proven that cryodeposits can form and remain stable on 
surfaces exposed to a space environment (refs. 6 to 8) . Such deposits could alter the 
thermal radiative properties of surfaces and may be critically important to the thermal 
heat balance of a spacecraft. 

These were the primary motivations for many of the more recent studies of the ra- 
diative properties of cryodeposits (refs. 9 to 15). Most of these investigations have been 
concerned with HgO and CC >2 cryodeposits at liquid-nitrogen temperatures and lower. At 
these low temperatures sublimation rates are negligible and lifetimes are effectively in- 
finite. Not only is the sublimation power loss from the surface of an HgO or a CC >2 cry- 
odeposit extremely small at 77 K, but it is several orders of magnitude smaller than the 
thermal blackbody radiation loss from a surface at that temperature. McConnell 
(ref. 15) considered the effect of sublimation on a cryodeposit that had formed on a 77 K 
substrate but did not define the effect that this term had on the overall heat balance. 

Little attention has been paid to the behavior of cryodeposits at higher temperatures 
where the effect of sublimation on the heat balance could be significant. As a result of 
higher temperatures, cryodeposits will sublime more rapidly and exhibit characteris- 
tics greatly different from those found at low temperatures. The particular concern 
here is how the overall equilibrium heat balance is affected when the sublimation heat 
loss becomes significant. Specifically, it is desirable to know when sublimation heat 
transfer is an important part of the overall energy balance and how long the cryodeposit 
will persist. 

In space flight these higher temperature cryodeposits could result from the jetti- 
soning of spacecraft life support products, the efflux of thrustor propellant, or the 
warming of substrates upon which the more common low-temperature cryodeposits ha/e 
formed. Also, if the presence of a cryodeposit would prove advantageous, means of 
growing the desired deposit could be developed. 

It is the purpose of this report to present a definitive study of the effect of sublima- 
tion on the heat balance of various cryodeposits and to consider the interaction between 
a subliming cryodeposit and the electromagnetic radiation to which it is subjected. In 
particular, a method of determining the surface temperature of the cryodeposit will be 
presented that will permit the evaluation of the relative importance of the various heat- 
transfe* mechanisms and the overall energy balance of an irradiated cryodeposit. Ex- 
pressions for the rate of mass loss from the surface of the cryodeposit permit good es- 
timates of the cryodeposit lifetime. As an example, numerical results are presented for 
one typical ice (H^O) cryodeposit. 
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SYMBOLS 

constant 

constant 

specific heat 

cryodeposit thickness 

time rate of change of deposit thickness 

free energy 

energy flux reaching substrate 
energy flux within the deposit 
energy flux incident on deposit surface 
energy flux entering the deposit 
thermal conductivity 
latent heat of sublimation 
molecular weight 
vapor pressure 

radiant energy absorbed per unit volume per time 

energy absorbed by the deposit 

net energy loss at the deposit surface 

radiation energy loss 

sublimation energy loss 

energy absorbed by the substrate 

gas constant 

reflectance of the deposit complex 

temperature 

surface temperature 

substrate temperature 

time 

sublimation rate 
absorbtance of the deposit 




e' effective emissivity ratio of the deposit 

p density 

cr Stefan Boltzmann constant 

r relaxation time 

Tj deposit lifetime 

ANALYSIS 
General Case 

Consider a planar cryodeposit that has formed on a thin metallic substrate in a 
vacuum. The deposit is subliming and there is a mass loss from the surface. The 
analysis presented herein starts with the existance of a cryodeposit and does not con- 
sider any further deposition. Actually, simultaneous deposition and sublimation (or 
evaporation) will exist to some extent in any cryodeposit. When deposition vapor pres- 
sure and substrate temperatures are such that the deposition rate exceeds the sublima- 
tion rate, a cryodeposit will continue to grow until an equilibrium deposit thickness is 
reached. At this time the deposition rate is equal to the sublimation rate. For the 
equilibrium situation to change requires that one of the governing parameters change. 
This could come about by a decrease, or even an effective elimination, of the deposi- 
tion vapor pressure or an increase in the substrate and hence the cryodeposit temper- 
ature . 

The net heat transfer at the surface of the cryodeposit is the algebraic sum of the 
sublimation energy loss , the thermal radiation energy loss , and the energy absorbed 
from any incident thermal radiation. The relative magnitude of these various heat- 
transfer mechanisms determines when they are importan , in the heat balance and when 
they can be neglected. 

The power loss, due to sublimation, from the surface of a solid in a high vacuum is 
a strong function of and is uniquely determined by the surface temperature. The details 
of the method used to determine this sublimation loss are given in the appendix. 

If the cryodeposit is opaque, or even if it is translucent but can be considered a 
"gray" emitter, the radiation loss Q rac j can be obtained simply and is written as 
e'er T^ where the effective emittance of the deposit has been denoted by e’ . 

The assumption of a gray deposit may not be very realistic for deposits that are 
highly transparent over large portions of the spectrum of the emitted radiation. How- 
ever, more exact solutions taking the spectral properties of the deposit into account re- 
quire detailed knowledge of these optical properties. And even with such detailed inf or - 
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mation, the actual computations could be very time consuming. McMahon (ref. 16), 
Gardon (refs. 17 to 19), and Merriam (ref. 12) treat this problem of the emissivity of 
translucent materials when the translucency and magnitude of this thermal radiation 
term for any particular cryodeposit are sufficient to warrant this refinement of the 
analysis. 

All the heat-transfer mechanisms are interrelated since each will affect the sur- 
face temperature and both the sublimation and radiation power losses are strong func- 
tions of the surface temperature. 

The analytical determination of the surface temperature of a subliming cryodeposit 
will now be considered with the additional assumption that the substrate is a constant 
temperature energy source. Furthermore, because of the temperature range under 
consideration and the small temperature variation through the deposit, it is assumed 
that there is no net radiant transfer between different regions of the cryodeposit. 

The temperature profile of a subliming cryodeposit is determined by the equation 
for the energy balance of a volume element within the cryodeposit: 


9T 


d 2 T 


pC n — = k ^ + q^ h 
p at .9 ^ Ds 


dx 


(i) 


o 

where p is the density (g/cm ), is the specific heat (J/(g)(K)), T is temperature 
(K), t is time (sec), k is the thermal conductivity (W/(cm)(K)), q a ^ s is the radiant 
power absorbed (W/cm 3 ), and the coordinate system is as follows: 



0-x 


and the boundary conditions are 


T(d) = T 


at 


~ ^sub^o^ + ^rad^o^ 


lx=0 


( 2 ) 

(3) 
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where Q gub is the sublimation power loss and Q ra( j is the radiation power loss from 
the surface of the cryodeposit. 

Equation (1) is a time -dependent equation that is quite difficult to solve because it 
contains a moving boundary, namely, the subliming deposit’s surface. Using a time- 
independent approach is justified if the fractional change in thickness of the solid over 
a typical thermal relaxation time is very small. Mathematically, this is expressed as 

r « 1 (4) 

o 

where d is the time rate of change of the deposit thickness and r is the thermal re- 
laxation time of the deposit. If this condition is satisfied, then it is possible to calcu- 
late temperature distributions and sublimation rates at any given instant and subse- 
quently determine the deposit thickness as a function of time. 

Thermal relaxation times for problems of this type are typically given by refer- 
ence 20 



. 4d 2 C p 

T < p — £• 


(5) 


n 


In the present report, deposit thicknesses <1 centimeter are considered and d are 
obtained from the time independent solution. Inequality (4) was found to be satisfied for 
all cases considered, and in fact dr/d never exceeded 10 -2 even for the irradiated de- 
posits. Thus, the analysis is continued using the steady -state form of equation (1): 


d 2 T 

dx 2 


_( *abs 


( 6 ) 


with the boundary conditions (2) and (3) . 

The power absorbed by an elemental volume q abg may be determined as follows . 

If L nc is the incident power per unit area at the deposit surface, a fraction (1 - #)I inc 

is removed from the incident radiation field where 9t is the reflectance of the deposit 

substrate composite. This absorbed power is divided into a portion that is transmitted 

through the deposit and is absorbed directly by the substrate Q__ and a portion that is 

ss 

absorbed throughout the deposit thickness . Thus, 

(1 - ^inc = Q ss + Q dep ™ 
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and 

f* ^abs = Q dep (8) 

y o 

The calculation of 9t, Q gg , Q dep , and q abg is a problem in the optics of thin ab- 
sorbing films and is not a thermophysical problem. For specular surfaces and non- 
scattering deposits, q abg and Q dep may be obtained using the equations of classical 
thin film optics (ref. 14); for nonspecular, scattering, and absorbing media, the method 
of Merriam (ref. 12) may be used. 

For convenience set 


‘labs = Q dep f ( x > 


so that 



(9) 


( 10 ) 


and f(x) expresses the fractional distribution of radiation absorbed within the deposit. 
We now solve equation (6) , which is rewritten as 

k nr ’ <W (x) (11) 

dx z 

Integrating equation (11) twice and using the boundary conditions (eq. (2)) and 
(eq. (3)) yield 


T = T - Q . 
o s ^sub 


(T.) H. - Q 
° k 


rad 


(T J - + 


Q, 


dep 


/V 

•'0 J 0 


f(x’)dx' 


( 12 ) 


This is the equation that determines the surface temperature of the cryodeposit. 
Equation (12) expresses the fact that the surface temperature of the deposit is lowered 
by the sublimation and radiation energy loss and raised by the energy absorbed from 
the incident radiation. Whether T is greater than or less than T depends on the 
relative magnitude of these three terms. Since the expression for Q gub (given in the 
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appendix) contains T q in both an exponential and a square root term, equation (12) is 
transcendental. The last term of equation (12) considers the distribution of the absorbed 
energy. This term is independent of T and T insofar as the optical constants of the 
deposit and the substrate are independent of temperature. Thus, the energy absorbed in 
the deposit is coupled to the thermophysical problem only through the fact that the de- 
posit thickness d changes due to sublimation. Any f(x) that expresses the fractional 
distribution of the absorbed energy, whether from a single pass through the deposit, 
from multiple reflections within the deposit, or even from a spectral analysis of the 
incident radiation, can be used in the solution of equation (12) if it can be specified. 

Since the absorbed power term is independent of temperature, equation (12) is alge- 
braic in nature and can be solved numerically. The absolute and the relative magnitude 
of the terms of equation (12) will determine which of the heat -transfer mechanisms has 
a negligible effect and which should be included to determine surface temperature. 

The analysis thus far applies if substrate temperatures are independent of time. 

This is a valid assumption for deposits that have formed on cryogen tanks. However, 
for other surfaces such as shadow shields, a deposit might play a major role in deter- 
mining substrate temperature . The energy loss due to sublimation will tend to lower 
the surface temperature. 

However, the assumption of constant substrate temperature may not be too restric- 
tive for the following reason: the thermal response time of the deposit -substrate com- 
plex will be of the order of, at least, seconds, while the characteristic thermal relaxa- 
tion time of the deposit (eq. (4)) will usually be in the millisecond range. Thus, one 
may use the time -independent solution to calculate an instantaneous surface temperature 
and power loss for any given substrate temperature and then use these quantities in a 
numerical solution of the equations that determine the temperature of the substrate. 

That is to say , the deposit layer will equilibrate fast enough to follow the more slowly 
changing substrate temperature. 


Optically Thick Case 

A fairly simple example has been chosen to illustrate the application of these 
methods. 

First, we assume that the incident radiation is absorbed by a translucent cryo- 
deposit according to the usual exponential law of absorption with no multiple reflections 
within the deposit. Thus, the intensity of the energy flux at any location within the de- 
posit is given by 


!f( x ) = I Q e 


■ax 


(13) 
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where I Q is the intensity of the entering incident energy flux at the surface and the total 
power absorbed in a single pass through the deposit is 


«dep - ' h - V 1 - e '“ d > 


(14) 


Since the power absorbed in an incremental thickness is 


= I ae _orx 
Ax ° 


(15) 


the fractional energy distribution within the deposit which is f(x) is given by 


d(I f ) 

. = f(x) =- “ 


-Q!X 


I o- I a 


1 - e 


-ad 


(16) 


Then, 


/ */ 


f(x')dx’ = 


1 - e 


-ad 


1 + JL ( e ~ ad - 1) 

ad 


(17) 


Thus, the absorbed power term in equation (12) for radiation absorbed according to the 
usual exponential law of absorption with no reflection at the substrate will be 


I. £ 1 + -i (e‘“ d - 1) 

° k ad 


(18) 


And for the optically thick case where the value of ad is large enough so that essen- 
tially no radiation reaches the substrate, this expression reduces to 


I 


o 



( 19 ) 


In the limit, the magnitude of the absorbed power term will be I d/k which is the 
opaque case. 


Since an opaque layer yields the maximum effect on the surface temperature, this 
limiting case will be used in the remainder of this example . In the case of the ice (H^O) 
frosts considered in reference 14, any thickness greater than 20 micrometers is op- 
tically thick. For such deposits 1^ = 0 since no radiation reaches the substrate and 

(1 " ^hnc = Q dep- Since (1 " (ref - 14) ’ Q dep = X inc’ and (eq ' (12)) becomes 


T o = T s-«sub< T o)j-^ T oJ +I incJ 


( 20 ) 


Equation (20) must be solved for T Q to obtain the net power lost by the irradiated de- 
posit 


Q *(T ) = Q , (T ) + e'uT 4 - I. (21) 

met' o' sub' o' o inc ' ' 

which may be positive or negative depending on the magnitude of the incident radiation 
and T g . 

Finally, the lifetime of a subliming deposit of thickness d is given by 

/ d 

d" 1 d(d) (22) 

o 

where the rate of decrease of the thickness d may be obtained from equations (Al) 
and (A 2) as 


l-'S.- gab (T o> (23) 

P P^sub 

RESULTS 

The sublimation power loss for water, carbon dioxide, argon, krypton, xenon, and 
oxygen over the temperature range that these materials are solids has been calculated 
and the results are shown in figure 1 . These calculations were based on equations (Al) 
and (A2) . The thermal radiation loss from the deposit is also shown in figure 1 for de- 
posit effective emittances of 1.0 and 0.1. Such a range should encompass the emittance 
of almost any cryodeposit and therefore yield an overall, qualitative comparison be- 
tween sublimation and radiation heat transfer. Figure 1 also indicates the level of the 
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Figure 1. - Sublimation loss and radiation loss as functions of surface temperature for water, 
(H 2 0), carbon dioxide (C0 2 l, xenon (Xe), krypton (Kr>, argon (A), and oxygen (0 2 ). 


power gained from a 300 K blackbody source and from 1 solar constant if all the inci- 
dent radiation is absorbed at the surface. 

Figure 1 shows that the heat balance analysis of a cryodeposit complex can be di- 
vided into high, middle (or transition), and low temperature regions. In the high- 
temperature region the deposit's thermal blackbody radiation energy loss is insignifi- 
cant compared with the sublimation energy loss . In the low temperature region the 
sublimation energy loss is negligible compared with the thermal blackbody radiation 
energy loss . Only in a fairly narrow temperature range are the sublimation and the 
thermal blackbody radiation heat losses of comparable magnitudes so that both must be 
considered in the heat balance analysis. When incident radiation is absorbed by the de- 
posit, it is possible for both the sublimation and the thermal blackbody power loss to 
become unimportant in a heat balance, particularly at low temperatures. 

Equation (20) has been numerically evaluated for T Q by the Newton-Raphson tech- 
nique for ice layers from 0.001 to 1 centimeter thick on substrates for which T = 190, 

s 

o 

200, 210, . . . 270 K for one particular set of frost parameters (p = 0.917 g/cm ; 

0^= 0.7 J/(g)(K); k / 0.04 W/(cm)(K); and e' = 1.0), both with and without an incident 
energy heat load from a 300 K blackbody. The effect of ice layers of thickness <10 mi- 
crometers may be obtained by simple extrapolation of the curves to be discussed. How- 
ever, it should be recognized that the assumption that the layer is optically thick (no 
radiation reaches the substrate) may not be valid for smaller thicknesses. The radia- 
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tionless case may be extrapolated without error. 

The temperature drop across the subliming frost layer is presented in figure 2 as 
a function of frost thickness for all substrate temperatures. For this particular ice 
deposit there is essentially no temperature gradient through the frost layer for a sub- 
strate temperature of less than 200 K. However, above this temperature the surface 
temperature of a 1 -centimeter -thick frost is appreciably lower than the substrate tem- 
perature, in fact, as much as 45° lower for the 270 K substrate case. 

The radiation equivalent to that coming from a 300 K blackbody tends to lower the 
temperature drop across the frost, but only in the case of the 190 K substrate is the 
effect of this irradiation enough to counteract the losses due to sublimation and radia- 
tion. The maximum effect that an incident radiation can have on the surface tempera- 
ture of a cryodeposit will occur at low temperatures where only the last term of equa- 
tion (20) must be considered. Thus the maximum temperature rise through a 
1 -centimeter deposit will only be 1 K for the 300 K blackbody case and 3. 5 K for 
1 solar constant absorbed at the deposit surface. If the deposit is not a blackbody ra- 
diator, the temperature drop will be diminished even further, but the effect is small. 

For the frosts considered here, eliminating the thermal radiation term entirely changes 

the surface temperature by only 0.18 K for the most sensitive case (T = 190 K, 

s 

d = 1 cm) . 

With the surface temperature of the frost known, it is now possible to solve equa- 
tions (22) and (23) for the frost lifetimes that have been plotted as a function of the frost 
thickness (fig. 3). Note that the curves for the lowest temperature substrates are quite 
straight and have slopes close to 45°. This is because there is little temperature gra- 
dient across these deposits so that the surface temperature is nearly independent of frost 
thickness and nearly equal to the substrate temperature. Consequently, the sublimation 
rates are nearly constant. As the temperature of the substrate increases, the upward 
curvature of the lifetime curve increases. This curvature for the frosts with an appre- 
ciable temperature gradient is due to the fact that thicker frost implies lower surface 
temperature, lower sublimation rate, and hence longer lifetime. If the temperature 
gradient through the frost had not been considered for the 270 K substrate temperature 
deposit, the 800-second life predicted would be only 17 seconds. 

The lifetime curves for the irradiated frosts are seen to be only slightly below the 
curves for the unirradiated frost, indeed, for frost temperatures >230 K the irradiated 
case is not indicated because they coincide, within plotting accuracy, with the unirra- 
diated lifetime curves . This behavior results because the radiant energy heat load of 
a 300 K blackbody is very small compared with the sublimation energy loss, which com- 
pletely dominates the energy transfer in this temperature range. 

Finally, note that the thickness as a function of time may be obtained from figure 3 
by starting time at some value of and moving down to the left on the appropriate 
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Net power loss from cryodeposit complex, 




curve. For instance, let us say that at t = 0, d = 0. 1 centimeter and T g = 230 K. 

Then = 110 seconds corresponds to t = 0, and we find that after 80 seconds 
(110 - 30) the frost thickness has decreased to 0.033 centimeter. 

The net power loss (eq. (20)) from the frost is plotted in figure 4 for the same range 
of d and T that was used in figure 2. The sublimation power loss decreases for in- 

fa 

creasing d since the surface temperature decreases. For T g > 230 K the radiant 
heat load makes a negligible contribution to the energy balance. The irradiated frost 
continues to lose power until T < 200 K. For T = 190 K the irradiated frost gains a 
net 0.026 watt per square centimeter, which is essentially independent of thickness. 

For T c « 190 K, the sublimation and radiation power losses are negligible, and the 
frost complex keeps the full 0.046 watt per square centimeter that it absorbs from the 
incident radiation field. 

If the radiation power loss is eliminated from the problem, by setting emittance e’ 
equal to zero, the change in the net power is even less than the effect of the incident 
energy heat load, for the frost and temperatures considered here, and is not significant 
for the higher temperature frosts. 


CONCLUSIONS 

The sublimation energy loss of a cryodeposit is quite sensitive to the surface tem- 
perature of the deposit. There is a narrow band of temperatures over which the sub- 
limation and radiation power losses are of the same magnitude, and both must be con- 
sidered in an energy balance. Above these temperatures sublimation becomes the 
dominant energy loss mechanism, and the thermal blackbody radiation loss may be ig- 
nored. In fact, under many conditions even an incident thermal (300 K) heat load be- 
comes insignificant compared with the sublimation energy loss. Consequently, the 
"absorptance" of the deposit complex is of little concern because the absorptance of in- 
cident radiation makes little difference in the power balance. At temperatures below 
this band, the sublimation energy loss is negligible and only radiation energy loss need 
be considered in a heat balance. 

Although the radiation power loss could be important in the low temperature region, 
results indicate that this loss has no appreciable effect on the surface temperature of a 
cryodeposit. Hence including the radiation energy loss in the surface temperature de- 
termination, is an unncessary refinement. 

Cryodeposits, and in particular ice deposits, that form on cold surfaces in a space 
environment, will persist for long times. Thermal, or even solar radiation should 
prove ineffective in removing a cryodeposit rapidly. In fact, for the higher substrate 
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temperatures, the incident energy load equivalent to radiation from a 300 K blackbody 
has essentially no effect on the calculated lifetime of the frost. Hence, it appears that 
in space the only means of eliminating an undesired cryodeposit in a reasonable time, 
is to raise the substrate temperature. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 11, 1971, 

124-09. 
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APPENDIX - ENERGY AND MASS LOSS FROM A SUBLIMING CRYODEPOSIT 


Consider a planar cryodeposit with surface temperature T q in a vacuum. The de- 
posit is subliming and there is an energy and mass loss from the surface. The sublima- 
tion mass loss rate is found from the Langmuir equation (ref. 21): 


W(T q ) = 


P ( T o> 

17.14 



(Al) 


where P(T ) is the vapor pressure (in torr) of the subliming solid of molecular weight 
M and absolute temperature T (K) . This equation is generally restricted to values of 
P(T ) of 1 torr or less. 

The sublimation power loss is then 


< W T o> = W < T c> L sub 


W 

2 

(cm )(sec) 


(A 2) 


where L gub is the latent heat of sublimation (in J/g). 

Evaluation of these losses requires knowledge of the vapor pressure at the cryo- 
deposit temperature. 

The most reliable method for correlating experimental observations on vapor pres- 
sure data is based on the calculation of the free energy. Expressions for the free energy 
are given by Kelley (ref. 22) for most elements and many inorganic compounds. Denot- 
ing the free energy by AF°, tie vapor pressure in torr may be derived by means of the 
relation 


A P’O 

In P = - ==— + In 760 (A3) 

RT 

Vapor pressure data are also available in references 23 and 24. 

Kelley (ref. 22) also presents values for the latent heat of sublimation in calories 
per mole which are needed in equation (A2) . 

When expressions for the vapor pressure and latent heat of sublimation are unavail- 
able for any particular cryodeposit, a good approximation for these values may be ob- 
tained in the following manner. 

For vapor pressures below about 1 torr, it has been observed that for most materi- 
als. the vapor pressure can be represented as a function of T by the relation 
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where A and B are constants. 

This is certainly valid for a narrow range of temperature variation. Assuming that 
the vapor pressure of ice obeys such an expression between the values given by Kelley 
(ref. 22) for the temperature range of 170 to 220 K, the data are fitted by the expression 

In P = 23.907 - (A 5) 

T 


The heat of sublimation can then be obtained from Clapryon’s equation 


L 


sub R 


d In P 



(A6) 


where the slope d In P/d(l/T) can be calculated from equation (A5). 
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